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Abstract. 3:2 and/or 3:1 twin high frequency quasi-periodic oscillations (HFQPOs) have been measured for the three micro- 
quasars GRO J1655-40, XTE J1550-564 and GRS 1915-1-105. For a test particle orbiting a rotating black on a stable circular 
orbit there exist two different orbits at which the vertical and radial epicyclic oscillations are in either a 3:1 or 3:2 parametric 
resonance for any choice of the black hole angular momentum a. If the two orbits are required to be frequency commensurable 
Keplerian orbits there is only one solution for the two orbit radii and a. This model predicts that the microquasars have the same 
a, and it predicts their black hole masses on the basis of the measured HFQPOs in agreement with the dynamically determined 
masses. Application of this model to the Galactic Center black hole Sgr A* using the recently measured QPOs (Genzel et al., 
2003, Aschenbach et al., 2004) leads to a black hole mass of (3.28 ± 0.13) x 10® Mq, and the same a as for the microquasars. 
The possibility that all four sources have a — 0.99616 suggests that this value is the upper limit of a imposed by general 
relativity. The same value for the lower orbit radius and the same value for a are also suggested by an analysis of the general 
relativistic expression for the radial gradient of the orbital velocity, which changes sign in a narrow annular region around the 
lower orbit when a > 0.9953. 

Key words, accretion discs, QPOs, relativity, - Galaxy: center - Galaxy: nucleus - Galaxy: fundamental parameters - X-rays: 
individuals: Sgr A* - X-rays: individuals: GRO J1655-40 - X-rays: individuals: XTE J1550-564 ~ X-rays: individuals: GRS 
1915-1-105 -X-rays: individuals: XTE J 1859-226 -X-rays: general 



1. Introduction 

High-frequency qasi-periodic oscillations (HFQPO) have been 
observed in quite a number of X-ray binaries containing a 
compact object, i.e. neutron stars or black holes. HFQPOs 
are considered to originate in the accretion disk around the 
compact object. A special group is obviously set up by the 
galactic microquasars GRO J1655-40, XTE J1550-564 and 
GRS 1915H-105. They are special because they do not show 
just one frequency but twin frequencies at a fixed ratio of 
3:2 (Abramowicz & Kluzniak, 2001, Remillard et al., 2002), 
which are 450, 300 Hz (Strohmayer 2001a, Remillard et al., 
1999), 276, 184 Hz (Remillard et al., 2002) and 168, 113 Hz 
(Remillard et al., 2003, McClintock & Remillard, 2004), re- 
spectively. Out of those three objects XTE J 1550-564 might 
have an additional power spectral density peak at 92 Hz, the 
significance of which is modest, however (Remillard et al., 
2002). 

Around the time the discoveries were made, it was sug- 
gested independently that some resonance effect in the accre- 
tion disk could give rise to such a 3:2 frequency behaviour 
(Kluzniak & Abramowicz, 2001a, b). There are three funda- 
mental oscillation modes for a test particle orbiting a black 
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hole, which produce the orbital or azimuthal (Kepler) fre- 
quency, the radial (in the orbital plane) epicyclic frequency 
and the vertical or polar (perpendicular to the orbital plane) 
epicyclic frequency. Initially it has been suggested that the 
3:2 frequency ratio is due to a coupling of orbital and ra- 
dial frequency (Abramowicz & Kluzniak, 2001, Schnittman 
& Bertschinger, 2004). However, the proponents of the res- 
onance theory then concentrated on a parametric 3:2 res- 
onance between the polar and the radial epicyclic modes 
(Kluzniak & Abramowicz 2002, Abramowicz et al. 2002, 2003, 
Abramowicz & Kluzniak 2003, Kluzniak & Abramowicz 2003, 
Rebusco, 2004), and Remillard et al. (2002) found solutions for 
a 3:2 frequency ratio for each of the three possible couplings, 
i.e. orbital/radial, orbital/polar and radial/polar. The method of 
associating a frequency to one specific oscillation mode has 
in principle the benefit that from frequency measurements the 
mass M, the angular momentum a (the KeiT parameter) of the 
black hole and the radial position of the orbiting test particle 
in the accretion disk can be determined. But with the measure- 
ment of even several frequencies the relevant equations remain 
underdetermined and more information is needed. So, in gen- 
eral the mass, which has been determined from the dynamics 
of the binary orbit for each of the three microquasars, has been 
made use of and a value of a could be determined, which, how- 
ever, is different for each of the three coupling possibilities. 
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Furthermore, as Remillard et al. (2002) have pointed out on 
their findings on XTE J 1550-564 that the frequency ratio might 
not be just 3:2 but that 3:1 is another possibiUty, one ends up 
with a choice for a of one out of six, still adopting M from dy- 
namical measurements. There may be a possibility to resolve 
this ambiguity by introducing another constraint in the reso- 
nance model. 



I have investigated whether both a 3:2 and a 3: 1 parametric 
resonance between the vertical and radial epicyclic frequen- 
cies would exist, of course at two different orbits, but that these 
two orbits are commensurable orbits. The commensurability is 
meant in the traditional sense that the Kepler frequencies of the 
two orbits come in ratios of natural numbers. Actually, there 
exists such a configuration and there is only one solution, so 
that the Kepler frequency of the inner orbit is three times the 
Kepler frequency of the outer orbit. Furthermore, this solu- 
tion allows only one value for the angular momentum, which 
is a = 0.99616. The mass of the black hole AI is uniquely 
determined after a choice for one of the two possible orbits 
has been made. It is likely that the discriminator between the 
orbits is the mass accretion rate. In any case the two possi- 
ble values for AI differ by a factor of 1.5 only. The masses 
of the three black holes in GRO J1655-40, XTE J1550-564 
and GRS I9I5H-I05, predicted by this model on the basis of 
their measured HFQPOs, agree in each case with the dynami- 
cally determined masses within their measurement uncertainty 
range. I have applied this model also to the Galactic Center 
black hole, Sgr A*, making use of the one quasi-period recently 
published by Genzel et al. (2003) and the quasi-periods pub- 
lished by Aschenbach et al. (2004). Again, the predicted black 
hole mass agrees very well with the dynamically determined 
mass (section 3). The possible relevance of the 3:2 resonance 
model to the periodicities reported for Sgr A* has been men- 
tioned in the papers of Kluzniak, Abramowicz & Lee (2003) 
and Abramowicz et al. (2004). 



I have checked the general relativistic expressions for the 
motion of a test particle orbiting a rotating black hole using 
the canonical Boyer-Lindquist functions for some unexpected 
behaviour near the two orbits and in the angular momentum 
range around a = 0.99616. Energy and angular momentum 
of the test particle remain unconspicuous, but the orbital ve- 
locity does not. Whereas for any a < 0.9953 the orbital ve- 
locity is monotonically increasing with decreasing radius over 
the full range of the radius down to the marginally stable orbit, 
for a > 0.9953 the orbital velocity starts to decrease with de- 
creasing radius over a small radial range before it again rises 
with decreasing radius. This annular range of decreasing or- 
bital velocity happens to be located where the 3:1 parametric 
resonance of the vertical and radial epicyclic frequencies ap- 
pears for a = 0.99616. To me it seems that this finding in not a 
chance coincidence but is backing the resonance model above. 
The detailed physics still have to be explored. Further details 
are given in section 4. Some discussion about the implications 
is presented in section 5. 



2. The model 

There are three fundamental cyclic gravitational modes asso- 



ciated with black hole accretion disks (Nowak & Lehr lI998l 
Merloni et al. 1999), which are the Kepler frequency (Ok), 
the disk perturbation frequencies in vertical and radial direc- 
tion called vertical (fly) and radial (fl^) epicyclic frequency. 
Each frequency depends on the central, black hole mass A'l, its 
angular momentum a and the radial distance r from the center. 
Equations^tolSlshow the relations, for which the standard no- 
tation of c=G=l is used. Physical length scales are in units of 
GM/c-^ and angular frequencies il are in units of c'^/GM. r = 
1 is defined as the gravitational radius rg. The conversion from 
observed linear frequency i' measured in Hz to the normalized 
il is obtained through i/ = rixc^/(2 n GM) with c, G, M in cgs 
units. 



S2y = Sij^ (1 — ■^pjj2 + ^3") 



nl^niii-l + ^-s^) 



(1) 

(2) 
(3) 



The model I propose to explain the twin HFQPOs observed 
at a ratio of 3:2 or 3:1 should have both fly/il^ = 3:2 and 
rty/flYi =3:1 for the same angular momentum a. Of course, 
this requires two different orbits with different radii called r32 
and Tai, which follow from Equations ^ to Figure [T] shows 
r32 and r^i versus (1 — a). There is no limitation on a, and 
r32 and r^i exist for any a. The second requirement is that 
the two orbits should be commensurable orbits, which in tra- 
ditional sense means that flKir = r^i) = n x i^K{r ~ '"32) 
in it's simplest form, with n a natural number. Figure |2] shows 
ilK{r = r3i,a)/nK{r = 7-32,0) versus (1 — a). The maxi- 
mum value of this ratio is 3.26 and the minimum value is 2.026. 
Therefore there is a solution but only one solution for n = 3, 
which in turn fixes the values of a to a — 0.99616, r^i = 1.546 
and r32 = 3.919. Values of a, r^i and r32 are truncated. 

Tabled summarizes the values of the relevant frequencies 
fi including the Lense-Thirring precession frequency f^LT = 
r^K — ^v (Merloni et al., 1999) for r3i and r32. The parametric 
resonances are for ratios of 3:1 and 3:2. I note that if a beat 
frequency of f^b = ^v — Or exists, which may be indicated 
by the observation of the 92 Hz frequency in XTE J 1550-564, 
then the parametric resonances would be a 3:2:1 triplet both at 
r^i and r32. Therefore I include Ob in Tabled 

It is evident from Table [J that fly is the upper frequency 
both in the 3:2 and the 3:1 resonance between fly and Or ne- 
glecting Ok. Accordingly, the black hole mass can be calcu- 
lated for the 3:1 resonance (M31) and the 3:2 resonance {AI32) 
from Equations0]and|5lbased on the HFQPO measured upper 
frequency Vup in Hz: 



M31/M0 = 4603.3/i^„p 
A/32/M0 = 3046.2/i^„p 



(4) 
(5) 



Equations 0] and |5] illustrate that AI31 and M32 differ by 
factor of 1.51 only. 
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Fig. 1. Radii of orbits for which flv/i^B. = 

riy/riij = 3:2 (dashed line) versus (1 ~ a). 
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Fig. 2. Ratio of Kepler frequencies J7k(?' — '"31, a)/i^K('' = 
'"32, a) versus (1 — a). 

Apart from the parametric resonances enforced by the ini- 
tial requirements, i.e. ilY{r — r3i)/i7R(r ~ r^i) = 3:1, 
^v{r = r32)f^R{r = rga) = 3:2 and flKir = r3i)/f7K(r = 
'"32) = 3:1 it turns out that there are 'cross' resonances con- 
necting rai and rs2, which are, apart from the Kepler com- 
mensurable orbits, ilji{r = ?'32)/^r(?' — ?'3i) = 4:3 and 
^Y^r = r32)/SlR(r = r^i) = 2:1 (c.f. Table^. In summary, 
^31 and r32 show an extraordinary high degree of parametric 
resonances. 

With a ~ 0.99616 some radii, useful to know, can be cal- 
culated. The event horizon is at rn ~ 1.088; the static limit 
in the equatorial plane is at rg = 2 and the marginally stable 
orbit is at Vms = 1.25 (defined by flu ~ 0; for the defini- 
tion by Bardeen et al. (1972) rms — 1.306), which in either 



Table 1. Compilation of frequency values for a = 0.99616, 
rai = 1.546 and r32 = 3.919. 



Q. 



rai 



J-32 



Qk 0.342684 0.114228 

fiv 0.142370 0.0942112 

^R 0.0475116 0.0628118 

^b 0.0948583 0.0313993 

f^LT 0.200313 0.0200171 



case means that both rai and r32 are stable orbits, and that r^i 
is within the ergosphere whereas r32 is outside of it. In fact, 
''31 = {th + ''o)/2, which is mid-way between the event hori- 
zon and the static limit in the equatorial plane. 



3. Comparison with observations 

The model so far has made use of just frequency ratios not so 
much of the absolute frequencies. The test of the model would 
consist of a comparison with otherwise measured masses or 
angular momenta. For the three microquasars GRO J1655-40, 
XTE J 1550-564 and GRS 1915H-105 measurements of their 
masses have been made but there is so far no way to deter- 
mine the angular momenta, so the test is necessarily restricted 
to their masses. Masses are usually derived from the determi- 
nation of the orbit parameters in case of a binary. This holds 
for the three microquasars but also for the black hole Sgr A* 
in the Galactic Center which is orbited by the S2 (SO-2) star 
(Schodel et al., 2002, Ghez et al., 2003). Genzel et al. (2003) 
have reported a QPO period of 16.8 min observed in two near- 
infrared flares. Aschenbach et al. (2004) have claimed addi- 
tionally quasi-periods around 219 s, 700 s, 1150 s and 2250 
s, with the 1 150 s period being consistent with the NIR period. 
This set of quasi-periods was found in the power density spec- 
tra of one X-ray flare observed with Chandra (Baganoff et al., 
2001) and a second X-ray flare with XMM-Newton (Porquet 
et al., 2003). Interesting in this context is that the frequencies 
corresponding to the latter three quasi-periods are close to a 
3:2:1 ratio. Enforcing such a ratio sequence in a best fit, it 
turns out that such a ratio is consistent with the measurements. 
Accordingly, I include the Galactic Center black hole in this 
test and predict the mass of Sgr A*. 

3.1. GRO J 165 5 -40 

The frequencies in question are 450 and 300 Hz (Strohmayer 
2001a, Remiflard et al. 1999) with fiy = 450 Hz. According to 
Equation |5] the predicted black hole mass is Mbh = 6.76 ± 
O.li\/0. The relative uncertainty of the mass is the same as 
that of the HFQPO measurements of « 1.5%. Dynamical mass 
measurements have been reported by Orosz & Bailyn (1997) 
with Mbh — 7.02 ± O.22M0 and more recently by Greene et 
al. (2001) who obtained a mass range of Afs// — 5.8 — 6.8Mq. 
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The agreement between model prediction and observation is 
quite satisfactory. 

3.2. XTE J 1550-564 

Frequencies of 276 and 184 Hz (Remillard et al., 2002) have 
been measured. With Q,v = 276 Hz Equation|5]predicts a black 
hole mass of Mbh = 11-04 ± O.2A/0. The relative uncer- 
tainty of the mass is again due to the accuracy of the HFQPO 
measurements of « 2%. Dynamical mass measurements have 
been reported by Orosz et al. (2002a) with a ±lcr mass range 
of Mbh = 8.4 — 11.6Mq, which is nicely matched by the 
predicted mass. 

3.3. GRS 1915+105 

The relevant 3:2 pair at 168 and 1 13 Hz has only recently been 
found (Remillai-d et al., 2003, McClintock & Remillard, 2004). 
With J7y = 168 Hz Equation |5l predicts a black hole mass of 
Mbh = 18. 13±O.36M0. The relative uncertainty of the mass 
is due the accuracy of the HFQPO measurements of « 2%. 
Dynamical mass measurements have been reported by Greiner 
at al. (2001) with Mbh = 14.0 ± 4.0A/o and Harlaftis & 
Greiner (2004) with Mbh = 14.0 ± 4.4Af0. Also in this 
case the predicted mass matches the dynamically mass within 
the ±1(7 range, but the measurements of this source illustrate 
nicely the potential uncertainty in dynamical mass measure- 
ments. Harlaftis & Greiner (2004) point out that the biggest 
uncertainty is related to the inclination of the orbit, which for 
GRS 1915H-105 is taken from the inclination of the associated 
jet assuming that the jet axis is orthogonal to the orbital plane 
of the binary. If, for instance, the inclination angle is changed 
from the adopted 66 to 56 the estimated black hole mass is 
Mbh = 16.9 ± 5.9Af0 (Harlaftis & Greiner, 2004). Actually, 
Kaiser et al. (2004) suggest that the inclination angle is 53 
putting the best estimate for the dynamical black hole mass at 
17.8 Af0, which is within the ilc measurement of Harlaftis & 
Greiner (2004), but it is also pretty close to the predicted value 
of 18.1M0. 

3.4. XTE J 1859-226 

Cui et al. (2000) report the measurement of a single HFQPO 
at 187 Hz associated with XTE J1859-226. Filippenko & 
Chornock (2001) announced a mass function of the binary of 
7.4Af0, and Brocksopp et al. (2002) discovered associated ra- 
dio jets suggesting that XTE J 1859-226 is another galactic mi- 
croquasar Orosz et al. (2002b) and McClintock & Remillard 
(2004) finally quote an estimated mass of the black hole be- 
tween 7.6 and 12.0 A/0. So far there is just one HFQPO mea- 
sured so that a prediction based on the proposed model is 
difficult, but there are only two options. If the 187 Hz fre- 
quency corresponds to the vertical epicyclic 3:2 resonance the 
predicted mass (Equation O is about 16.3 Af0 and a second 
HFQPO at 125 Hz, i.e. the radial epicyclic frequency, should be 
looked for If, on the contrary, the 187 Hz frequency is the result 
of the radial epicyclic frequency the expected vertical epicyclic 



frequency is at 281 Hz and the predicted mass is about 10.9 
A/0, which fits quite nicely in the mass range given by Orosz 
et al. (2002b) and McClintock & Remillard (2004), but these 
authors mark the mass range as quite uncertain. A deep search 
for a second HFQPO would be extremely useful. 



3.5. The Galactic Center black hole SgrA* 

Aschenbach et al. (2004) have grouped the frequencies indi- 
cated in the power density spectra of two NIR flares (Genzel et 
al. 2003) and two X-ray flares. There are three groups for which 
there is evidence that each of the frequencies of one group oc- 
curred in at least one NIR and both X-ray flares. These three 
groups of quasi-periods in seconds, ordered by NIR, Chandra, 
XMM-Newton, are group 1: (733, 692, 701), group 2: (1026, 
1117, 1173) and group 3: (-, 2307, 2178); the duration of the 
NIR flare was too short to cover the longest period. Looking 
at the average quasi-period within one group it appears as if 
the ratio of the average quasi-periods follows a 1:1.5:3 relation 
which is equivalent to a 3:2:1 relation in frequencies. This is at 
least consistent in view of the frequency measurement uncer- 
tainty. Each of these eight quasi-periods has been observed at 
low wave number which means that the systematic uncertainty 
of each period determined from the associated power density 
spectra is up to ±20%. Enforcing a 3:2:1 ratio sequence in a 
best fit the highest frequency is i^up = 1.4035 ± 0.0574 mHz 
or a period of 712±29 s. The two remaining quasi-periods are 
at 1069±44 s and 2138±87 s. In the paper of Aschenbach et 
al. (2004) we associated the quasi-period of group 1 with the 
vertical epicyclic oscillation at the marginally stable orbit and 
the quasi-period of group 2 with the radial epicyclic oscillation 
at the radius at which the radial oscillations have their maxi- 
mum frequency. This together with the association of the 219 s 
quasi-period to the Kepler period of the marginally stable orbit 
led to a unique value for the mass and the angular momentum 
of the black hole. For the quasi-period of group 3 (2150 s) we 
had no convincing suggestion. This together with the fact that 
we needed to employ different radial positions for the different 
frequencies is unsatisfactory. This is taken care now by the as- 
sociation of the quasi-periods of group 1 , group 2 and group 3 
with the epicyclic modes at just one radius. 

Unlike the microquasars which show two frequencies 
SgrA* appears to show three frequencies of about similar 
power spectral densities (Aschenbach et al., 2004), so that it 
makes sense to classify this as an example of the 3:1 reso- 
nance for which Equation @] applies, and the black hole mass 
of SgrA* is predicted to Mbh = (3.28 ± O.13)xlO^M0. 
Measurements of the orbital motion of the star SO-2 (S2) 
around the GC black hole have resulted in Mbh = 3. 7 ±1.5 
X IO6M0 (Schodel et al., ImT^, Mbh = 4.07±0.62 x 
1O^M0 (Ro/8kpc)3 (Qjjg^ et al.. 2003 1 with Rq the distance to 
the b lack ho le, and Mbh = 3.59±0.59 x IO^Mq (Eisenhauer 
et al.,"2004"l. Errors quoted for each of the three mass measure- 
ments are la. The mass predicted by the present model agrees 
with both the Schodel et al. (2002) mass and the mass given by 
Eisenhauer et al. (2003) but it is slightly out of the la range of 
the Ghez et al. (2003) measurement for R = 8kpc. 
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Application of Equation |5] checking for a 3:2 resonance, 
gives Mbh = (2.17 ± 0.1)xlO'^Mo, which is ruled out by 
the dynamical measurements. This supports that in the frame- 
work of this model the accretion disk around Sgr A* is in a 
3:1 resonance. The third frequency associated with the reso- 
nance state 2 may be either a beat frequency, i.e. ily-^R or the 
first harmonic of 57r. In either case a fourth frequency could 
be expected with is either fiy+^R or the third harmonic of 
J7r. This would correspond to resonance state 4, and the ex- 
pected period is 534±22 s. Such a period seems to be present 
in the power density spectrum of the Chandra flare at 494±41 
s and in the NIR flare of June 15, 2003 at 498±100 s. It is ab- 
sent in the NIR flai-e of June 16, 2003 and the XMM-Newton 
flare (Aschenbach et al., 2004). With a and Mbh given the 
orbital (Kepler) period is Pxir = r^i) = 296 ± 12 s; at 
the marginally stable orbit Pk^t — rms) — 242 ± 10 s. 
The Chandra flare shows increased power spectral density at 
^Chandra = 256 ± 10 s and the XMM-Newton flare has signif- 
icant power at Pxmm = 219 ± 15 s (Aschenbach et al., 2004). 
These data are consistent with the view that the shortest quasi- 
periods observed in the X-ray flares correspond to Keplerian 
motion along orbits between r^i and r„is, the actual radius 
might be different from flare to flare, not so are Q.y and Q,ji, 
of course. 

4. A physical mechanism? 

With the distinct values for a, r^i and r^2 given it comes to 
mind to look for a physical reason why these quantities are spe- 
cial. I checked the general relativistic expressions of energy, an- 
gular momentum and orbital frequency for a test particle mov- 
ing around a rotating black hole on a stable orbit. These ex- 
pressions (Shapiro & Teukolsky, 1983) change monotonically 
with r in (a, r) space. But this is not the case for the orbital 
velocity v^"^^ described in the ZAMO-frame (the Zero Angular 
Momentum Observer or Bardeen observer), v^^' is represented 
by Equation|6l(Muller & Camenzind, 2004): 



t;(*) ^^in-w). 



(6) 



The Boyer-Lindquist functions are (e.g., Miiller & 
Camenzind, 2004): 



P\/A 



A = r^ - 2Mr + a^ 



p^ = r'^ + a^cos^i 



£2 = (^2 + a'^f - a'^Asin'^0, 



(7) 

(8) 

(9) 

(10) 

(11) 

(12) 



9 is the polar angle with 6 = tt/2 defining the equatorial 
plane. 

As Miiller and Camenzind (2004) state one usually assumes 
Keplerian rotation for r > Vms and sets $7 = Q,k (Equation 



u 


— 


— stnO 
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Fig. 3. Orbital velocity w*^*^ versus orbital radius for a — 
0.99616. Outside the radial range shown w'^*^ decreases mono- 
tonically with radius. 

[0. Setting 9 = tt/2 and il = Hk I have solved numerically 
Equation|6lover the full (a, r) space. For large values of r and 
small values of a, i;'*^ behaves as expected from Newtonian 
mechanics, i.e. w'*^ increases monotonically with decreasing 
r. But this monotonic behaviour changes for a > 0.9953; w*-*^ 
develops a maximum at small values of r, decreases with de- 
creasing r until it reaches a minimum and from thereon in- 
creases again with decreasing r until r^s is reached. Figure|3] 
shows an example for a = 0.99616. The depression of w'*) of 
0.2% is small but it increases with increasing a. Figure|4]illus- 
trates how the radial position of the local maximum (rmax) and 
the local minimum (rmin) of v'*' change with a. The plot also 
demonstrates that always rmin > fms- Such a variation with r 
is absent for the particle energy E and its angular momentum 1. 
A main role, that this is different for the orbital velocity, seems 
to play the appearance of oj, the frame-dragging frequency or 
potential for angular momentum, in Equation|6] The change of 
t)(*^ implies that for Tmin < r < r^ax dE/dv'^'^^ < and 
dl/dv^'^^ > 0. 

According to standard notion the orbits between rmin and 
rmax are stable. The region rmm < r < rmax is the only re- 
gion with dv^'^> /dr > 0; it also has 9^w*-*^/9r^ = 0. It is ba- 
sically the rate of change of the orbital velocity per unit length 
of r seen by the ZAMO. It has the dimension of a frequency. I 
therefore define an angular 'frequency' ilc ~ 27r g^ \max, 
by which the radial position is defined at which dv^"^^ /dr 
has the maximum value for a given a. Figure |5] shows Vtc in 
comparison with Q,ji, Q,v and Q,k versus (1-a). It might be 
a chance coincidence but Vtc has the same value as Q.i^ for 
a — 0.99616 at r3i, indicated by the lowest cross in Figure|5l 
It appears that there is a 1:1 resonance between 51c and 0,^, 
and this is the only 'resonance' with Q.r, Vty and VLk for 
a < 0.998, which is the 'Thorne' limit (Thorne, 1974). In sum- 

a (■*) 



mary, Q.c = 27r- 



dr 



picks the same values for a and r as 
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Fig. 4. As a function of (1 — a) are shown the radial position 
Tmax of the local maximum of u^*) (curved dashed line), the 
radial position rmin of the local minimum of w^*) (curved dot- 
ted line). The vertical dashed line is at a = Oc = 0.9953. The 
cross marks rsi and a = a/ ~ 0.99616, which is the solution 
of the parametric resonance model with commensurable orbits. 
The solid line represents the innermost marginally stable orbit 
Tms defined by ilji = 0. 

the parametric resonance model with commensurable Kepler 
orbits. 



5. Discussion 

The model predicts the mass of each of the four black hole ob- 
jects with well estimated masses quite satisfactorily from mea- 
sured HFQPOs. Nevertheless, there are a couple of interesting 
coincidencies which I would like to point out and discuss. 

It is interesting to note that the model based on paramet- 
ric resonances between the vertical and radial epicyclic fre- 
quencies constrained by the requirement that any associated 
orbits should be commensurable orbits has only one solution, 
which selects just two orbits with radii r^i and r32, and restricts 
the angular momentum of the black hole to just one value of 
a = Qf — 0.99616. Even more interesting is the result that 
the orbital velocity w'*^ for the ZAMO changes its otherwise 
monotonical behaviour with radius and shows a local minimum 
and a local maximum which define a narrow radial annulus 
with dv^'^^ / dr > for a > 0.9953 = Oc. To me it seems that 
this feature of very rapidly rotating black holes has been over- 
looked so far. This certainly needs more study. Furthermore, 
the fact that ilc = '^'^^^bF' \max^ ^n{r ^ r^i), suggests that 
there is a link between the 'strange' behaviour of v^'^^ and the 
radial epicyclic oscillation, such that 51c(a/) = ^ii:(?'3i, a/) 
is the condition to excite the radial and vertical epicyclic os- 
cillations. I note that so far all this is produced in a single test 
particle environment. 



Fig. 5. As a function of (1 — a) are shown flKir — fmim a) 
and Vlxir = r„iax,ci) (label K), which makes up the upper 
parabola like contour; the contour labelled V is for riy(r = 
rmin, a) and ily(r — rmax,a) and the lowest parabola like 
contour (R) is for ftR{r = rmm^a) and i}R{r = rmax,a). The 
vertical dashed line is for a = af = 0.99616 and the crosses 
mark the corresponding frequencies found from the commen- 
surability considerations. Their positions are clearly within the 
radial section bounded by r„ 



belled c represents ftc = 2tt 



dr 



The model predicts that there is just one single value of 
a = Qf, and the HFQPOs of the three microquasars and the 
GC black hole are consistent with this. The same value of a 
for four sources suggests that this is some sort of final state 
of a; if correct this would replace the former 'Thorne' upper 
limit of a = 0.9982 - 0.9978 (Thorne, 1974). This would also 
imply that the excitation of the radial and vertical epicyclic os- 
cillations at rai somehow regulate the angular momentum flow 
towards the black hole to an equilibrium point, rai is located 
in the ergosphere halfway between the horizon and the static 
limit; the distance Sr to the ergosphere limit is Sr — 0.45 in the 
equatorial plane. A radial oscillation with an amplitude greater 
than Sr would traverse the ergosphere boundary back and forth. 

The oscillations at r^i are likely to be powered by mass 
accretion and one could envisage the following scenario. At 
low mass accretion rates both the vertical and radial oscilla- 
tions are stable and energy is stored in the corresponding wave; 
with increasing accretion rate the vertical oscillations become 
unstable because they are not any longer capable to store the 
energy and/or angular momentum but the energy is released 
in flares or bursts. I suggest that this is the state in which the 
GC black hole is. Observations show that the quiescent state 
X-ray luminosity of Sgr A* is about 1 1 orders of magnitude 
below the Eddington limit (Baganoff et al., 2001, Porquet et 
al., 2003), which might indicate that the accretion rate is fairly 
low, although other mechanisms for the low luminosity have 
been proposed. The time scale for the rise time of the flare 
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should be a half up to a few times of the period of the verti- 
cal oscillations, i.e. about 350 to 1000 s, and the duration of 
such a flare should be about a few times of the period of the 
radial oscillations, i.e. 40 to 120 min. These timescales have 
been observed in the Chandra (Baganoff et al., 2001) flare and 
the XMM-Newton flare (Porquet et al., 2003, Aschenbach et 
al., 2004). Recently Bower et al. (2004) have reported measure- 
ments by which Sgr A* has been spatially resolved at a wave- 
length of 7 mm, and they estimate itssize to be <2.4x 10^^ cm 
at a wavelength of 1.3 mm for a distance of 8 kpc. I note that 
the diameter of the region, where the vertical epicyclic mode 
is at work, is 2xr3i = 1.5 x 10^^ cm, whereas the diameter of 
the r32 annulus is 3.8x10^^ cm, which appears to be incom- 
patible with the wavelength extrapolated radio measurements, 
if the X-ray emitting region is not more extended than the radio 
source. This result seems to confirm that the Sgr A* emission 
region is restricted to an area with a radius r < 7-31. 

At even greater accretion rates there are no stable vertical 
oscillations any more, the 3:1 resonance disappears, and proba- 
bly a significant fraction of the energy of the accreting matter is 
lost in vertical direction in a non-periodic fashion or stored in 
the radial mode, the amplitude of which will grow. Eventually 
the radial wave will reach r32 and power the vertical and radial 
modes at rs2 because of their parametric 'cross' resonances 
(section 2). This is the 3:2 state I suggest the three microquasars 
are in. The X-ray luminosities of the three microquasars are 
very high, close to or possibly slightly above the Eddington 
limit. The accretion rate is certainly fairly high. Eventually, out- 
bursts may also occur in the 3:2 state, when the vertical mode 
oscillations become unstable temporarily. Of course, a quan- 
titative relation between oscillation state and accretion rate is 
still to be worked out. 

The three microquasars are known to have jets, and Bower 
et al. (2004) state that their radio measurements of Sgr A* are 
consistent with jet models, so that all four sources might have 
jets. Therefore one might ask the question, whether the prop- 
erties twin or triple HFQPOs, a ^ af and existence of a jet 
come together all three, and vice versa. For instance, has an 
accreting black hole with a jet a = a/ and HFQPOs? If so, 
it could be profitable to search for instance the extragalactic 
jet sources, make an estimate of the HFQPOs from Equation 
131 and Equation [S] and when found determine the black hole 
mass to percentage accuracy. Finally, though the idea might 
be foolish, it comes to mind. Is the instability of the vertical 
epicyclic oscillation at r = r^i the seed for the creating of ac- 
cretion driven black hole jets? The decrease of w^*) of about 
2% between r^ax and rmin would correspond to an apparent 
loss of orbital kinetic energy of about 4.3 MeV/nucleon in non- 
relativistic Newtonian approximation. Whether and how this 
energy can be made use of is unclear. 

6. Conclusion 

I have constructed a model built upon parametric 3:1 and 3:2 
resonances between the vertical and radial epicyclic oscilla- 
tions of a test particle orbiting a rotating black hole. If these 
two orbits are required to be commensurable orbits as far as 
their Keplerian frequencies are concerned there exists only one 



solution, which specifies uniquely the orbit radius r^i = 1.546 
of the 3 : 1 resonance, the orbit radius r32 = 3 . 9 1 9 of the 3 :2 res- 
onance and the angular momentum a/ = 0.99616. The analy- 
sis of the general relativistic expression for the orbital velocity 
v*^*^ in the ZAMO frame shows a deviation from the mono- 
tonic relationship between u^*^ and r for a > 0.9953 = Uc- 
dv^"^^ /dr changes its sign over a narrow radial annulus, and 



there Q.r — 211- 



c= ri/j. It seems that this condition is 



responsible for the excitation of radial and vertical epicyclic os- 
cillations. The effect of dv^"^^ /dr changing its sign is likely to 
be related to energy and angular momentum transfer via verti- 
cal and radial waves. This effect seems to have been overlooked 
in the past and needs more study. Since so far the model is 
based on test particle motion the entire effect is due to general 
relativity and for rapidly rotating black holes. It is not a col- 
lective effect. This, in some sense, consistent model predicts a 
new upper limit of a — af — 0.99616, which is slightly lower 
than the 'Thorne' limit (Thorne, 1974). The model predicts the 
mass of fully spun-up black holes on the basis of measured 
twin (3:2) or triple (3:1) HFQPOs to the accuracy of the fre- 
quency measurements. It is suggested that the 3:1 resonance is 
active at low accretion rates, and the 3:2 resonance can be ob- 
served at very high accretion rates. At intermediate accretion 
rates Q,ii{r — rsi) may be the only persistent frequency. The 
comparison of the masses derived with this model, on the basis 
of the measured HFQPOs, with the masses derived from dy- 
namical measurements, i.e. binary orbit characteristics, shows 
an excellent agreement for the three microquasars GRO J 1655- 
40, XTE J 1550-564 and GRS 1915H-105 and lowers the mass 
uncertainty significantly. The frequencies measured recently by 
Aschenbach et al. (2004) for Sgr A* appear to be arranged in a 
3:2:1 ratio and the derived mass, which is 5.7 orders of magni- 
tude higher than the lowest mass of the microquasars involved, 
agrees quite well with the dynamical mass at a much reduced 
uncertainty. 
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